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Regarding a d—dimensional spherically symmetric line element in the context of Einstein-A grav¬ 
ity, the hydrostatic equilibrium equation of stars is obtained. Then, by using the lowest order 
constrained variational (LOCV) method with the AV18 potential and employing microscopic many 
body calculations in the modern equation of state, the structure properties of neutron stars are in¬ 
vestigated. Regardless of cosmological point of view and considering arbitrary positive and negative 
values of the cosmological constant, the maximum mass of the neutron stars and their corresponding 
radius in 4-dimensions are computed. The results show that there is an upper limit for the maximum 
mass of neutron star for positive cosmological constant (Mmax < 1.68Mq). On the other hand, it 
is shown that the Einstein gravity cannot explain the structure of neutron star with negative A. 

Other properties of neutron stars such as; the Schwarzschild radius, average density, compactness 
and Buchdahl- Bondi bound are studied. In addition, by using the Buchdahl-Bondi bound for neu¬ 
tron stars, stability of these stars is investigated. Finally, the dynamical stability is investigated 
and shown that the neutron stars follow the dynamical stability in this gravity. 


I. INTRODUCTION 


Most of the solar system phenomena such as the Mercury precession can be successfully explained in the context of 
Einstein gravity, but when we intend to study beyond the solar system or strong gravity regimes, this theory encounters 
some problems. Accelerated expansion of our universe is one of the main challenging subjects of cosmologists EH, 
whereas Einstein gravity cannot properly explain it. In order to interpret this expansion, some various candidates 
have been proposed by many authors. Among various alternatives, one may point out the modified gravities like 
Lovelock gravity , brane world cosmology ^-[l^ , scalar-tensor theories Elnl3 ^-nd F{R) gravity [2ll - l3ll| . 

From the cosmological point of view, the dark energy is an unknown form of energy which is hypothesized to 
permeate the spacetime, tending to accelerate the expansion of the universe [32| . The dark energy is the most accepted 
hypothesis to explain the observations since the 1990s indicating that the universe is expanding at an accelerating rate. 
Another model of these interesting candidates is related to the cosmological constant. In this theory, the cosmological 
constant may add to the Einstein Lagrangian [s^ . [s^ for explaining such acceleration. The cosmological constant 
is the simplest possible form of dark energy, which leads to the current standard model of cosmology known as the 
A-CDM model and provides a good fit to many cosmological observations. 

From the astrophysical point of view, stars reach the equilibrium state due to the balance between gravitational 
force and internal pressure. Therefore, in order to study the structure of stars, we should obtain a suitable hydrostatic 
equilibrium equation (HEE) and then we can study the structure properties of stars. The neutron and quark stars 
are in the category of celestial objects with high mass, high density and small radius, and therefore, they are called 
compact objects. Due to this fact, we need to take into account the effects of general relativity (spacetime curvature) 
for investigating their structures. The first HEE of stars in 4-dimensional Einstein gravity has been studied by Tolman, 
Oppenheimer and Volkoff fTOVi J35l-K^ . Also, the physical characteristics of stars using TOV equation have been 
investigated by many authors [38l - l45l| . Recently, survey of compact objects in different gravity has attracted the 
attention of researchers, in particular, the models of modified gravity which involve higher curvature terms or a scalar 
field [dbl - ls^ . On the other hand, if one is interested to study the structure and evolution of stars in different theories 
of gravity with various sources, the new HEE must be obtained. Therefore, in recent years, the generalizations and 
modifications of TOV equation have been investigated in various theories of gravity, such as: dilaton gravity [s^, 
gravity’s rainbow [s^, massive gravity [H^, F{R) and F{G) gravities [dll ? and p ost-Newtonian theory (see 
also Refs. [ddl - fT^ for more details). 

The outline of our paper is as follows. In the next section, we consider a spherical symmetric metric and obtain the 
HEE in the Einstein-A gravity for arbitrary dimensions (d > 3). Then, we investigate stability and energy conditions 
for the equation of state of neutron star matter in 4-dimensions. Also, we study effects the cosmological constant on 
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maximum mass and its radius of the neutron stars. Next, we examine effects of the cosmological constant on other 
properties of neutron stars such as; Schwarzschild radius, average density, compactness, Buchdahl-Bondi bound and 
dynamical stability. We finish our paper with some closing remarks. 


II. HEE IN THE PRESENCE OF COSMOLOGICAL CONSTANT 

As first step, we regard d—dimensional Einstein—A gravity to obtain its related HEE. The action of the Einstein 
gravity with the cosmological constant in arbitrary dimensions, d, is given by 

Ig = -^ [ d^xV^{R-2A) + lMatt, ( 1 ) 

Jm 

where R and A are, respectively, the Ricci scalar and the cosmological constant in d—dimensions and iMatt is the 
action of matter field. Varying the action m with respect to the metric tensor 5 ^, the equation of motion can be 
written as 

R; + \Rg; + Ag; = KdT;, ( 2 ) 

where Gd is the gravitational constant in d-dimensions and define as Gd = GVd- 4 , where G denotes 

the four dimensional gravitational constant and Vd -4 is the volume of extra dimensions. Also, the constant c is the 
speed of light in the vacuum. In addition, R'^ and are the symmetric Ricci tensor and energy-momentum tensors, 
respectively. 

Here, we want to obtain the static solutions of Eq. ©. For this purpose, we assume a spherical symmetric spacetime 
in the following form 


ds^ = f{r)df - ^ - r'^dfll, (3) 

where /(r) and g(r) are the unknown metric functions of radial coordinates and dfl^ is the line element of unit 
(d — 2 )-dimensional sphere 


d- 2 i-l 

dnl = del + '^Y[ dOf ■ (4) 

i=2 j=l 

On the other hand, the energy-momentum tensor for a perfect fluid is 

= {c^e + P)U^U'^ - Pg^'', (5) 

where e and P are density and pressure of the fluid which are measured by local observer, respectively, and is the 
fluid four-velocity. Using Eqs. © and ([ 5 ]) and the metric introduced in Eq. we can obtain the components of 
energy-momentum for d-dimensions as follows 

To° = ec^ & Tl = T| = Ti = ... = = -P- (6) 

We consider the metric m and Eq. (jS]) for the perfect fluid and obtain the components of Eq. ([5]) in the following 
forms 


zv 2 2 A 2 , ^2d3 , . d2 , 

KdC r e = Ar + [I - g) - —rg , 


( 7 ) 


Kdfr^P = -Ar^f-^^{l-g)f+^rgf, (8) 

AKdfrP = -iArf - ^^ 3^4 (1 - g) f ^ ^ ^ ^ 25 /"] /, (9) 

r 

where di = d — i and also, /, g, e and P are functions of r. We note that the prime and double prime are, respectively, 
the first and second derivatives with respect to r. 
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Using Eqs. 0 - ® and after some calculations, we obtain 

dP {ch+p)r 

dr ^ 2/ 


Now, we obtain /' from Eq. as follows 


2 [r^ (A + KdP) + ^ (1 - g)] / 

rgd2 

Then, we calculate g{r) by using Eq. 0 with the following form 

, , , 2A 2 c^KdMT (f) 

dld2 (i27'‘'^2/2j.d-3 ’ 


( 10 ) 


( 11 ) 


( 12 ) 


where M = J 1" 1® gamma function, in which r(l/2) = r(l) = 1 and r(a; + 1) = xr(a;). 

By considering Eqs. (HU and (ED, and inserting them in Eq. ED, we can extract the HEE of Einstein-A gravity 
in d-dimensions as 


dP 

dr 


^ (A + %KaP) 




r( 




2 ir 


3572 


c^KdM - ^r<^^ 


j{ch + P), 


(13) 


where for 4-dimensional case (d = 4), Eq. (1131) reduces to the following equation 

dP [3c^GM + r3(Ac^ + 12^GP)] _ 

dr c^r [6GM — c^r {Ar^ -I- 3)] ^ 

Also, Eq. (fT4l) leads to primitive TOV equation for vanishing A [ssl - fs^ . It is evident that for d = 3, Eq. (fT^ 
reduces to the HEE obtained by Diaz Isfl. As a result, we have obtained a global form of HEE for arbitrary dimensions 
(d > 3) in Eq. (ED- 


III. STRUCTURE PROPERTIES OF NEUTRON STAR 
A. Equation of state of neutron star matter 

Obtaining equation of state for neutron stars enables one to study its properties. The interior region of the neutron 
star is a mixture of neutrons, protons, electrons and muons. These constituents are in charge neutrality and beta 
equilibrium conditions (beta-stable matter) [sij . Recently, the microscopic constrained variational calculations based 
on the cluster expansion of the energy functional has been employed to obtain the equation of state of neutron star 
matter |^, [s^- The two-nucleon potentials which are used in these studies are the modern Argonne AVIS [s^ 
and charged dependent Reid-93 [s^. One of the advantages of these methods is that there is no need for any free 
parameter in formalism, and they result into a good convergence. Also, the results in this method are more accurate 
comparing to other semi-empirical parabolic approximation methods. This is due to a microscopic computation of 
asymmetry energy which is carried on for the asymmetric nuclear matter calculations. In fact, it is necessary to have 
a microscopic calculation with the modern nucleon-nucleon potentials which is isospin projection {T^) dependant [s^ . 
Here, we are interested in using the lowest order constrained variational (LOCV) method with the AV18 potential 
[821183| for calculating the modern equation of state for neutron star matter and investigating some physical properties 
of neutron star structure. We use the LOCV method to calculate the energy of our system. The LOCV method is a 
useful tool for determination of the properties of neutron, nuclear and asymmetric nuclear matter at finite and zero 
temperature [STl - f^ . It is notable that, the LOCV method is a fully self-consistent formalism and it does not bring 
any free par ameters into calculation. It employs a normalization constraint to keep the higher order term as small as 
possible ( 91 L l92j| . The functional minimization procedure represents an enormous computational simplification over 
unconstrained methods that attempt to go beyond lowest order. 

A suitable trial many-body wave function is 


-0 = F(l), 


( 15 ) 
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FIG. 1: Equation of state of neutron star matter (pressure, P (10^® erg/cm®) versus density, e (lO'^® gr/cm®)). 

where (j) is the uncorrelated ground-state wave function of N independent neutrons, and is a proper A^-body 
correlation function. Applying Jastrow approximation [o^, F is rewritten by 

F = Sl[fitj), (16) 

i>j 

where 5" is a symmetrizing operator. Also, we consider a cluster expansion of the energy functional up to the two-body 
term 


Em = 


1 

N m) 


— El + E2, 


The energy per particle up to the two-body term is 


(17) 


Em = Ei+E2, (18) 

2 (i) ^ ^ j 

where Ei = X]i=-i- - | 2 m ^~r E 2 = ~ P) one-body and two-body energy term, respec¬ 

tively. The operator i^(12) has been given in Ref. (see Refs. [11,191,1^ for more details). Our obtained equation 
of state of neutron star matter is shown in Fig. [T] 

In order to obtain more details for equation of state of neutron star matter introduced in this paper, we study both 
stability and energy conditions in the following subsections. 


1. Stability 

In order to investigate the stability of the introduced equation of state for a physically acceptable model, one expects 
that, the speed of sound {v) be less than the speed of light (c) [^, [^. This stability condition can be written in the 
following explicit form 


0<u2= <c2. (19) 

By using the above condition and Fiq. [1] and comparing them with diagrams related to sound speed-radius 
relationship plotted in Ref. [9^ (see Fig. 2 in Ref. [9^, for more details), it is evident that this equation of state 
satisfies the inequality 0 < < c^. Therefore, the introduced equation of state is stable. 
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TABLE I: Energy conditions for neutron stars. 


ep {l(fkg/cm^) Pp (10^kg/cm^) NEC WEC SEC DEC 
3895.99 1910.39 / / / / 


2. Energy Conditions 

One of interesting properties of the mentioned equations of state is related to its central density. It is notable that 
the obtained central density of this equations of sta te is about 3.9 x 10^^ g cm~^ and this density is larger than the 
normal nuclear density ep = 2.7 x lO^^g cm~^ |lOCll |. 

Now, we want to study the energy conditions. In this regard, we investigate null energy condition (NEC), weak 
energy condition (WEC), strong energy condition (SEC) and dominant energy condition (DEC) at the center of the 
neutron stars. These energy conditions can be written in the following explicit forms 


NEC 

-)> 

-Pc + £c ^ 0, 


(20) 

WEC 


Pc + Ec P 0, 

& Ec > 0, 

(21) 

SEC 


Pc + Ec P 0, 

& 3Pc -|- Ec > 0, 

(22) 

DEC 

-)■ 

Cc > |Pc| , 


(23) 


where Pc and Cc are the pressure and density at the center of neutron star (r = 0), respectively. Considering Fig. [T] 
and the above conditions (l^ - (l^ . our results are presented in Table HI 

According to the Table U we observe that all the energy conditions are satisfied. Therefore, the equation of state 
of neutron star matter introduced in this paper are suitable and consistent with all versions of energy conditions. 

The equation of state of neutron star matter satisfied both stability and energy conditions. Now, we focus on 
investigation of maximum mass and radius for the neutron stars in Einstein-A gravity. 


B. Properties of neutron stars in Einstein-A gravity 


Considering the maximum gravitational mass of a neutron star for stability against collapse into a black hole, one is 
able to recognize differences between neutron stars and black holes. In other words, there is a critical maximum mass in 
which for smaller than such mass, degenerate pressure originated from the nucleons prevents an object from becoming 
a black hole [8l|. Therefore, obtaining the maximum gravitational mass of neutron stars is of great importance in 
astrophysics. Due to many errors in direct ways of measuring the neutron star mass by observations of the X-ray 
pulsars and X-ray b ursts , one is not able to obtain an accurate mass for these stars. On the other hand, using the 
binary radio pulsars fioMo^ . one may obtain highly accurate results for the mass of neutron stars. 

Now, by employing the equation of state of neutron star matter presented in Fig. [T] and numerical approach for 
integrating the HEE obtained in Eq. m, we obtain the maximum mass of neutron stars. The neutron star mass 
and radius depend on the central mass density (cc). We can consider the boundary conditions P(r = 0) = Pc and 
m(r = 0) = 0, and integrate Eq. (fT51) outward to a radius r = i? in which P vanishes for selecting an Ec- This leads 
to the neutron star radius R and mass M = m{R). 

In Ref. [9^, the Einstein gravity without cosmological constant has been investigated, and maximum mass of 
neutron stars has been obtained using the modern equations of state of neutron star matter derived from microscopic 
calculations. It was shown that the maximum mass of neutron star is about I.GSMq. In the present paper, we are 
going to consider the cosmological constant and obtain the maximum mass for neutron stars by employing the modern 
equations of state and study other properties of these stars. Our results are presented in Tables El and HTH 

The value of the cosmological constant and its fine tuning is open questions for scientists. In addition, its consistent 
values in various scenarios are different and therefore, it is allowed to regard it as a free parameter. On the other hand, 
from the standpoint of cosmology its order of magnitude may be approximately 10“®^ m“^. Such order of magnitude 
may be valid in a special cosmological model with large scale structure. Meanwhile, in order to regard the effects of 
A in local scales (such as near the black holes, neutron stars and other massive objects), it is not necessary to follow 
the fine tuning values of A is large scales. Here, we regard the value 10“^^ m~^ and also other values (at least as a 
toy model) to obtain the effects of cosmological constant on the neutron star structure. According to our results, the 
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TABLE II: Einstein-A gravity for 4-dimensions. 


A(m ^) 

(Mn) 

R(km) 

Rsch{km) 

e (10*® (? cm °) 

»7(10“*) 

z 

Mbb (Mq) 

0 

1.68 

8.42 

4.95 

1.34 

5.88 

0.56 

2.54 

1.00 X 10“®° 

1.68 

8.42 

4.95 

1.34 

5.88 

0.56 

2.54 

1.00 X 10“®^ 

1.68 

8.42 

4.95 

1.34 

5.88 

0.56 

2.54 

1.00 X 10“®'* 

1.68 

8.42 

4.95 

1.34 

5.88 

0.56 

2.54 


TABLE III: Einstein-A gravity for 4-dimensions. 


A(m ^) 

(Mr,) 

R{km) 

Rsch{km) 

e (10*®g cm ®) 

r,(10“*) 

z 

Mbb (Mq) 

1.00 X 10“** 

0.78 

6.65 

2.30 

1.26 

3.46 

0.24 

2.00 

5.00 X 10“*^ 

1.12 

7.47 

3.30 

1.27 

4.42 

0.34 

2.25 

1.00 X 10“*^ 

1.56 

8.25 

4.60 

1.32 

5.57 

0.50 

2.49 

5.00 X 10“*® 

1.62 

8.34 

4.77 

1.33 

5.72 

0.53 

2.51 

1.00 X 10“*® 

1.67 

8.40 

4.92 

1.34 

5.86 

0.55 

2.53 

1.00 X 10“** 

1.68 

8.42 

4.95 

1.34 

5.88 

0.56 

2.54 

1.00 X 10“*® 

1.68 

8.42 

4.95 

1.34 

5.88 

0.56 

2.54 


cosmological constant has no significant effect when we consider its value identical to 10“®^ (see Table HIl for 
more details). Therefore, in order to examine its effects, we should consider a (an ad hoc) model in which its value is 
less than 

Based on Table m considering the cosmological constant as a positive quantity, the maximum mass of neutron 
star decreases as A increases. In other words, the maximum mass of neutron stars with the positive values of 
cosmological constant are in the range Mmax < I. 68 M 0 . Also, when we consider the cosmological constant to be very 
small (A < 10“^^), the maximum mass reduces to results obtained in the Einstein gravity [^. In other words, for 
A < 10“^'* this constant does not affect the maximum mass and other properties. So, considering the obtained value 
for the cosmological constant from the cosmological perspective (A ~ this quantity does not affect the 

structure of neutron stars. 

We plot diagrams related to the maximum mass versus central mass density (radius) in Figs. [2]and[31 The variation 
of maximum mass versus the cosmological constant is also shown in Fig. [5] for A > 0. 

As we can see, mass of neutron star versus radius changes when the A increases. Indeed, one of interesting results 
of this paper is related to change in diagrams of mass versus radius, so that by increasing A, these diagrams change 
from ordinary shapes of the neutron stars to those of quark stars (see Fig. [3]). In other words, by increasing A, we 
encounter with a change of state (phase transition) of neutron star to the quasi quark star. On the other hand, we 
consider negative values of the cosmological constant and plot the density and pressure versus radius of neutron star. 
As one can see, the maximum and minimum values belong to the center and surface of neutron star, respectively (see 
Figs. OandE]), but the behavior of diagrams related to mass-radius relation of these stars are not logical (see Fig. [7]). 
In other words, Einstein gravity with the negative values of the cosmological constant are not a suitable theory for 
explaining the neutron stars with the mentioned equation of state. 

Our results show that, by decreasing the cosmological constant (A < 10“^^m“^), the maximum mass and also 
the radius of this star are not modified. In other words, the cosmological constant does not significantly affect the 
maximum and radius of neutron stars for A < 10“^^ m~^. As a result of this paper, when the value of the cosmological 
constant is about 10 “®^ m~^, this constant does not play a sensible role in the structure of neutron stars, but by 
considering larger values for it (about A > 10“^^ m~^), the maximum mass and its radius are modified (reduced). 

In the following, we investigate other properties of neutron stars by considering the cosmological constant such as; 
the Schwarzschild radius, average density, compactness, Buchdahl-Bondi bound and dynamical stability. 
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FIG. 2: Maximum mass of neutron star versus central mass density tc (10^®gr/cm®), for A = 1 x 10“^“^ (continuous line), 
A = 1 X 10“^® (dotted line), A = 5 x (dashed line) , A = 1 x (dashed-dotted line), A = 5 x 10“^^ (bold line) and 

A = 1 X 10“^^ (bold-dotted line). 


1. Schwarzschild Radius 


Now, we obtain the Schwarzschild radius for the obtained masses in Einstein-A gravity in the following form 


RSch 


(^3GM + 9G'2M2^ A^c 

1/3 

C 


Ac 

r / 


\ -|l/3’ 


3GM + +9G^M^ A^c 


(24) 


where Rsch is the Schwarzschild radius, G = 6.67 x 10“^^ Nw?/kg'^ and c = 3 x 10® m/s. It is notable that, by 
applying the cosmological constant to the Einstein gravity, the Schwarzschild r adiu s modified and also, considering 
A = 0 in Eq. (??), it reduces to the Schwarzschild radius in Einstein gravity as [H 


2GM 

^Sch — n 


(25) 


To find the Schwarzschild radius of neutron stars, we use of Eq. (??). The results show that increasing the 
cosmological constant and decreasing the maximum mass of neutron stars, the Schwarzschild radius decreases (see 
Table ITTll for more details). 


2. Average Density 

The average density of the neutron stars in 4—dimensions can be written as the following form 

__ 3M 
47ri?® ’ 

where the results for various the cosmological constant are presented in Table IIIII 


(26) 


3. Compactness 


The compactness of a spherical object may be defined as the ratio of the Schwarzschild radius to the radius of 
object as 


V = 


RSch 

R 


(27) 
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FIG. 3: Maximum mass of neutron star versus radius for A = 1 x 10 (continuous line), A = 6 x 10 
A = 1 X 10“^® (dotted line) and A = 1 x 10“^^ (dashed-dotted line). 
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FIG. 4: Maximum mass of neutron star versus cosmological constant (xlO ^^) in 4-dimensions. 
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FIG. 5: Density of neutron star versus radius. 



r(cm) 


FIG. 6: Pressure of neutron star versus radius. 


which may be indicated as the strength of gravity. We obtain the values of rj in the Einstein-A gravity and find that 
increasing A leads to decreasing the strength of gravity. 


4- The Gravitational Redshift 


Using the equation (USD, we obtain the gravitational redshift in 4-dimensions for this gravity as 

1 

- - 1 . 


1 - 


2GM A p2 


(28) 


As one can see, decreasing A leads to increasing the gravitational redshift (see Table mil. The maximum value of 
the gravitational redshift is about 0.56. According to Table HID one finds variation of A changes the gravitational 
redshift, drastically. This is a way for adjusting the values of A based on observational evidences. 
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FIG. 7: Maximum mass of neutron star versus central mass density Cc (10^®gr/cm®), (left panel) and radius (right panel) for 
A = —9 X 10“^^ (continuous line) and A = —5 x (dotted line). 


5. Buchdahl-Bondi Bound 

According to Buchdahl-Bondi bound |l06l - ll0^ . for a static spherically symmetric perfect fluid sphere, maximum 
permitted M — R ratio should be M < ||t i?. For the Einstein-A gravity this bound modified as [l09j| 

M<Mbb (29) 

where Mbb = — ^^R^ + A-\/l + 3AR‘^. Using the above condition, we investigate the stability condition 

of neutron stars in Einstein-A gravity. The calculations are presented in Table IIIII The results show that obtained 
maximum mass of neutron stars for all values obtained in Table imi are stable. 


6. Dynamical Stability 


The dynamical stability of th e stellar model against the infinitesimal radial adiabatic perturbation was introduced 
by Chandrase khar in R ef. [ng. Next, this stability condition was developed and applied to astrophysical cases by 
many authors [UMil. The adiabatic index ( 7 ) is defined as 


ec^ -I- P dP 
c^P de ' 


(30) 


It is notable that, in order to have the dynamical stability, 7 should be larger than | (7 > | = 1.33) everywhere 
within the isotropic stars. So, we plot a diagram related to 7 versus radius for different values of the cosmological 
constant in Fig. [H] As one can see, this stellar model is stable against the radial adiabatic infinitesimal perturbations. 

On the other hand, in order to present more investigations for the inside of neutron stars in Einstein-A gravity, 
we plot density and pressure versus radius of this star in Figs. l^andfTOl These figures show that density (pressure) 
values have a maximum and a minimum in the center and surface of neutron star, respectively, as expected. 


IV. CLOSING REMARKS 

In the first part of this paper, we have considered a spherical symmetric metric and extracted the hydrostatic 
equilibrium equation of stars for d—dimensional Einstein-A gravity (d > 3). Then, we have investigated the maximum 
mass of neutron stars by using the modern equations of state of neutron star matter derived from microscopic 
calculations. The results showed that the cosmological constant affects the maximum mass of neutron stars. In other 
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FIG. 8: Adiabatic index of neutron star versus radius. 



FIG. 9: Density of neutron star versus radius. 


words, considering the A > 0, the maximum mass for neutron stars decreases when A increases. In this case, maximum 
mass of neutron stars is in the range M^ax < 1.68Mq. 

One of the interesting results obtained in this paper is the behavior of diagrams of mass versus radius of the neutron 
stars. By increasing A, these diagrams showed a transition from the neutron star to the quasi quark star. On the other 
hand, considering A < 0, the behavior of diagrams related to mass versus radius (central mass density) of neutron 
stars is not logical, so the Einstein gravity with the negative values of cosmological constant can not explain the 
existence of neutron stars. In order to investigate the interior of neutron stars in Einstein-A gravity in more details, 
we have plotted density and pressure versus radius of this star. The results showed that the density and pressure were 
maximum and minimum of values at the central and surface of neutron stars, respectively. 

Our results showed that, when the approximate value of the cosmological constant was 10“®^ m~^, this constant did 
not play a sensitive role in the structure of neutron stars, but by taking larger values for it (about A > 10“^^ 
the maximum mass and its radius were reduced. Another interesting result of this paper is related to the effect 
of cosmological constant on the so-called gravity strength. We have shown that when the positive values of the 
cosmological constant increases, the strength of gravity decreases, and this effect leads to decreasing the maximum 
mass of neutron stars. 

Finally, we have studied the stability condition for these stars using Buchdahl-Bondi bound and found that the 
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FIG. 10: Pressure of neutron star versus radius. 


obtained maximum mass of neutron stars for all values obtained in Table Hill were stable. Then, considering adiabatic 
index, we showed that these stars were stable against the radial adiabatic infinitesimal perturbations. 
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